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(57) Abstract 



A system and method for adaptive interpolation of image data selects one of a plurality of different interpolation kernel functions 
based on one or more variables, such as a selected format for the visible presentation of the interpolated image data on an output imaging 
media, a selected type of the output imaging media for visible presentation of the interpolated image data, and an appearance characteristic, 
such as sharpness or smoothness, of the visible presentation of the interpolated image data. The system and method select an appropriate 
interpolation kernel function to improve an appearance of the visible presentation of the interpolated image data on the output imaging 
media. 
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SYSTEM AND METHOD FOR ADAPTIVE 
INTERPOLATION OF IMAGE DATA 

Fidd of the Invention 

The present invention relates to image processing techniques, and, more 
particularly, to interpolation techniques for expansion of digital image data for 
presentation on output media. 
Discussion of Related Art 

A digital laser imager forms a viable representation of an image on 
imaging media, such as photosensitive film, based on a set of image data 
representative of the image. The image data contains a plurality of digital image 
values acquired by an input imaging device. In a medical imaging environment, 
the input imaging device may comprise, for example, a diagnostic device, such as 
a magnetic resonance (MR), computed tomography (CT), conventional 
radiography (X-ray), or ultrasound device. Thus, the image data may represent an 
image of a physiological object obtained for diagnostic purposes. Each of the 
digital image values corresponds to one of a plurality of pixels in the original 
image, and represents an optical density associated with the respective pixel. The 
laser imager processes the digital image values to generate laser drive values. 
Each of the laser drive values represents an exposure level calculated to 
reproduce, on the imaging media, the optical density of a pixel in the original 
image. The laser drive values are used to modulate the intensity of a scanning 
laser, thereby exposing imaging media with a particular level of exposure. The 
laser imager subsequently develops the imaging media to form the viable 
representation of the original image. 

The number of pixels represented by the image data can be significantly 
less than the resolution of the imaging media. For example, laser imagers include 
scanning lasers capable of producing an output resolution of 300 pixels per inch 
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on the imaging media. Depending on the resolution of the input imaging device, 
however, the original image may be represented by image data comprising as little 
as a 256 by 256 array of pixels. The laser imager may provide different formats 
that enable a user to select the number of images to be formed. The number of 
5 images formed on the imaging media determines the size of each image. If a 

particular format requires that the laser imager form an image in a two inch (4.08 
cm) by two inch (4.08 cm) area on the imaging media, a 512 by 512 pixel set of 
image data must be expanded by at least twenty-five percent to match the 300 
pixels per inch resolution of the laser imager. With formats requiring larger sizes, 

10 it is necessary to expand the image data by much greater amounts. An eight inch 
(20.32 cm) by ten inch (25.4 cm) format for an entire sheet of imaging media, for 
example, would require eight-times expansion of the 512 by 512 pixel set of image 
data, resulting in a 4096 by 4096 pixel set of image data. In contrast, formats 
producing smaller sizes require that the image data be reduced. 

1 5 Expansion involves resampling of the originally sampled image data at a 

higher rate dependent on the amount of expansion required. Reduction involves 
resampling of the originally sampled image data at a lower rate dependent on the 
amount of reduction required. Resampling is carried out by performing 
interpolation on the original image data to effectively fit a curve to the original 

20 image data. For expansion, the curve is resampled at a higher rate than the 

original image data to create intermediate pixel points between those available 
from the originally sampled image data. For reduction, the curve is resampled at a 
lower rate than the original image data. The resampling process transforms the 
image data into an image having a higher spatial resolution for expansion and a 

25 lower spatial resolution for reduction. Unfortunately, the resampling rates for 
expansion can degrade image quality by loss of high frequency information, 
introducing perceived blurring and other artifacts into the expanded image data. 
To preserve image quality, it would be desirable to have an interpolation function 
that exhibits an ideal low-pass filter spectrum. The ideal low-pass filter would 

30 remove replicates of the frequency introduced by the resampling. The 
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implementation of an ideal low-pass filter is not feasible, however, because h 
would require the availability of many surrounding data points for each 
interpolated pixel. 

Several proposed interpolation functions seek to approximate the effect of 
5 an ideal low-pass filter. Such interpolation functions generally fall into two 

categories, commonly referred to as (1) "two-point" interpolation and (2) "four- 
point" interpolation. One type of "two-point" interpolation involves "pixel 
replication" interpolation, whereby sample points between two adjacent original 
pixels are simply assigned the intensity values of the closest adjacent original pixel. 

1 0 Another type of <4 two-point" interpolation involves "bilinear" interpolation, 

whereby sample points between two adjacent original pixels are assigned linearly 
weighted intensity values based on distance from the closest adjacent original 
pixel. Other "two-point" interpolation techniques are implemented by third order, 
fifth order, seventh order, etc., polynomial kernel functions. 

1 5 "Four-point" interpolation may involve the application of a cubic spline 

function. The four-point interpolation carried out by application of the cubic 
spline function involves much more computation than two-point interpolation. 
The cubic spline function has been found to produce a better spatial frequency 
response than two-point interpolation functions in most cases, providing a closer 

20 approximation of the ideal low-pass filter spectrum. For certain types of image 
data, however, a cubic spline interpolation function may produce a less than 
optimum interpolated image, relative to two-point interpolation functions such as 
bilinear interpolation or pixel replication. To improve interpolation results for 
image data representing alphanumeric information, for example, it may be 

25 necessary to use a different type of interpolation than that used for data 
representing pictorial information. 

Some existing interpolation techniques enable switching between different 
interpolation functions based on the type of image data to be interpolated. United 
States Patent No. 5,008,752 discloses an interpolator for applying two different 

30 types of interpolation functions depending on the type of image data present. This 



PCT/US95/13560 



technique switches between different interpolation functions when a difference 
between successive pixel values exceeds a predetermined threshold, indicating a 
transition from pictorial to text information. United States Patent No. 4,468,693 
discloses an interpolation apparatus that allows a user to manually switch between 
cubic spline interpolation, bilinear interpolation, and pixel replication interpolation 
according to the user's preference. United States Patent No. 5, 327,257 discloses 
an apparatus that applies more than one interpolation function toaset of image 
data to achieve a desired appearance. This apparatus creates both a sharp 
interpolated image and a smooth interpolated image by applying a cubic spline 
function to a set of image data, and then creating a third image by combining the 
sharp and smooth interpolated images with selected weighting coefficients to 
produce a sharp or smooth resultant image. 

The above patent disclosures generally address the effectiveness of 
different types of interpolation functions for different types of image data. 
However, the effect of other image processing variables has not been adequately 
addressed. For example, the effectiveness of a particular interpolation function 
may vary according to the degree of expansion required by the format of the 
image data on the imaging media. If a format requires a large image expansion 
factor for an image, it is, of course, necessary to create a greater number of 
intermediate pixels between adjacent original unexpanded pixels. As a result, 
there will be a higher chance of introducing imaging artifacts such as blurring. If 
the format selection calls for an smaller expansion rate, however, the blurring will 
be less perceptible. Hence, an interpolation function considered inappropriate for 
use with larger expansion rates, due to high incidence of blurring, may 
nevertheless be appropriate for use with smaller expansion rates. 

In addition, the response characteristics of different imaging media can 
greatly affect the visible results of the interpolation. It has been observed, for 
example, that dry silver film has a development speed that is significantly slower 
than that of conventional silver halide film, in both toe and shoulder speeds. This 
slower development speed can render interpolation functions previously 
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considered effective for silver halide films less effective for dry silver film. The 
slower development speed requires that the scanning laser transmit a higher 
exposing energy to the imaging media to achieve a desired maximum optical 
density. The higher exposing energy translates to a larger dynamic range for the 
scanning laser in order to expose every gradation on the imaging media. 
However, the available dynamic range of existing, commercially available scanning 
lasers, such as laser diodes, may be limited for required exposing wavelengths in 
the infra-red regions. As a result, the scanning laser may be incapable of 
producing the desired number of gradations. The slower development speed of 
the dry silver film, given the limitations of the scanning laser, can result in loss of 
high spatial frequencies. The visible manifestation of the loss of high spatial 
frequencies may be a loss of apparent edge sharpness producing an undesirable 
blurring in high contrast areas. The blurring is particularly apparent in those areas 
involving a transition from image data representing alphanumeric text at a 
minimum density to image data representing pictorial information at much higher 
densities approaching maximum density. Therefore, an interpolation function 
considered appropriate for silver halide film may produce inadequate results with a 
dry silver film. 

As a further variation, presentation of the image data on alternate output 
media such as a cathode ray tube (CRT) display monitor, having a different 
response characteristic than film, can also affect the visible results of interpolation 
functions. A CRT monitor medium has several problems when compared to a 
print film for image presentation. For example, a CRT monitor suffers from 
severe dynamic range limitations (maximum intensity to minimum intensity range), 
relative to film. The dynamic range of a CRT monitor limits perceptible image 
gradation presentation to about 100- 150 levels as compared to print films, which 
may provide 1 500-2000 levels of gradation. Thus, like dry silver film, CRT 
monitors may suffer blurring in high contrast areas. It is noted, however, that the 
blurring suffered by a CRT monitor may be much more pronounced than that on 
dry silver film. In addition, a CRT monitor produces viewing glare and geometric 
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distortion not present in films. Thus, the perception of fine edges and detail is 
much more difficult for CRT monitors than for print films. Due to the smaller 
number of gradation levels, and the glare and distortion in a CRT monitor, an 
interpolation function appropriate for a film may be inappropriate for a CRT 
S monitor. 

Existing interpolation techniques generally fail to address the effect of 
format and/or media variations on interpolation results. Consequently, image 
quality continues to be a concern for interpolated images formed with different 
formats on different imaging media. Accordingly, there is a need for an 
1 0 interpolation system capable of selecting interpolation functions appropriate for 
particular formats and/or different imaging media. 

Summary of the Invention 

In view of the limitations of existing interpolation techniques, the present 

15 invention is directed to a system and method for adaptive interpolation of image 
data. The system and method of the present invention are capable of selecting one 
of a plurality of different interpolation kernel lookup tables based on a format 
selected for the image data. The system and method of the present invention also 
are capable of selecting one of a plurality of different interpolation kernel lookup 

20 tables based on an imaging media selected for presentation of the image data. By 
addressing variations in format and media, the system and method of the present 
invention enable application of an appropriate interpolation kernel to the image 
data to improve viable results. As a further advantage, the system and method of 
the present invention are capable of selecting one of the interpolation kernel 

25 lookup tables to achieve a particular appearance characteristic based on a degree 
of sharpness or smoothness desired by a user. 

Additional features and advantages of the present invention will be set 
forth in part in the description that follows, and in part will be apparent from the 
description, or may be learned by practice of the present invention. The 

30 advantages of the present invention will be realized and attained by means 
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particularly pointed out in the written description and claims hereof as well as in 
the appended drawings. 

To achieve the foregoing advantages, as broadly embodied and described 
herein, the present invention provides, in a first embodiment, a system for 
5 performing interpolation on image data representative of a plurality of pixels 
within an image to produce interpolated image data representing interpolated 
pixels, the system comprising media selection means for selecting one of a 
plurality of imaging media for presentation of the image, format selection means 
for selecting one of a plurality of formats for presentation of the image on the 

10 imaging media, appearance characteristic selection means for selecting one of a 
plurality of appearance characteristics of the presentation of the image on the 
imaging media, a memory storing a plurality of different interpolation kernel 
lookup tables, each of the interpolation kernel lookup tables containing a plurality 
of interpolation coefficients, interpolation kernel selection means for selecting, 

1 S based on a combination of the format selected by the format selection means, the 
media selected by the media selection means, and the appearance characteristic 
selected by the appearance characteristic selection means, one of the interpolation 
kernel lookup tables stored in the memory, and interpolating means for applying 
the interpolation coefficients contained in the selected interpolation kernel lookup 

20 table to the image data to produce the interpolated image data. 

In a second embodiment, the present invention provides a system for 
performing interpolation on image data representative of a plurality of pixels 
within an image to produce interpolated image data, the system comprising format 
selection means for selecting one of a plurality of formats for presentation of the 

25 image on imaging media, a memory storing a plurality of different interpolation 
kernel lookup tables, each of the interpolation kernel lookup tables containing a 
plurality of interpolation coefficients, interpolation kernel selection means for 
selecting, based on the format selected by the format selection means, one of the 
interpolation kernel lookup tables stored in the memory, and interpolation means 



7 
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for applying the selected interpolation kernel lookup table to the image data to 
produce the interpolated image data. 

In a third embodiment, the present invention provides a system for 
performing interpolation on image data representative of a plurality of pixels 
5 within an image to produce interpolated image data, the system comprising media 
selection means for selecting one of a plurality of imaging media for presentation 
of the image, a memory storing a plurality of different interpolation kernel lookup 
tables, each of the interpolation kernel lookup tables containing a plurality of 
interpolation coefficients, interpolation kernel selection means for selecting, based 

10 on the media selected by the media selection means, one of the interpolation 

kernel lookup tables stored in the memory, and interpolation means for applying 
the selected interpolation kernel lookup table to the image data to produce the 
interpolated image data. 

In a fourth embodiment, the present invention provides a method for 

1 5 performing interpolation on image data representative of a plurality of pixels 
within an image to produce interpolated image data representing interpolated 
pixels, the method comprising the steps of selecting one of a plurality of imaging 
media for presentation of the image, selecting one of a plurality of formats for 
presentation of the image on the imaging media, selecting one of a plurality of 

20 appearance characteristics of the presentation of the image on the imaging media, 
selecting one of a plurality of different interpolation kernel lookup tables based on 
a combination of the selected media, the selected format, and the selected 
appearance characteristic, each of the interpolation kernel lookup tables containing 
a plurality of interpolation coefficients, and applying the interpolation coefficients 

25 contained in the selected interpolation kernel lookup table to the image data to 
produce the interpolated image data* 

In a fifth embodiment, the present invention provides a method for 
performing interpolation on image data representative of a plurality of pixels 
within an image to produce interpolated image data, the method comprising the 

30 steps of selecting one of a plurality of formats for presentation of the image on 
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imaging media, selecting one of a plurality of different interpolation kernel lookup 
tables based on the selected format, each of the interpolation kernel lookup tables 
containing a plurality of interpolation coefficients, and interpolation means for 
applying the selected interpolation kernel lookup table to the image data to 
produce the interpolated image data. 

In a sixth embodiment, the present invention provides a method for 
performing interpolation on image data representative of a plurality of pixels 
within an image to produce interpolated image data, the method comprising the 
steps of selecting one of a plurality of imaging media for presentation of the 
image, selecting one of a plurality of different interpolation kernel lookup tables 
based on the selected media, each of the interpolation kernel lookup tables 
containing a plurality of interpolation coefficients, and interpolation means for 
applying the selected interpolation kernel lookup table to the image data to 
produce the interpolated image data. 

It is to be understood that both the foregoing general description and the 
following detailed description are exemplary and explanatory only, and not 
restrictive of the present invention, as claimed. 

Brief Description of the Drawings 

The accompanying drawings are included to provide a further 
understanding of the present invention and are incorporated in and constitute a 
part of this specification. The drawings illustrate exemplary embodiments of the 
present invention and together with the description serve to explain the principles 
of the invention. 

Fig. 1 is a functional block diagram of a system for adaptive interpolation 
of image data, in accordance with the present invention; 

Fig. 2 is a functional block diagram of an interpolation processor for 
incorporation in the system of Fig. 1, in accordance with the present invention; 
and 



WO 96/16380 



PCT/US95/13560 



Fig. 3 is a flow diagram illustrating a method for adaptive interpolation of 
image data, in accordance with the present invention. 

Detailed Description of the Preferred Embodi ments 
S Fig. 1 is a functional block diagram of a system 10 for adaptive 

interpolation of digital image data, in accordance with the present invention. As 
shown in Fig. 1, system 10 includes a format selector 12, an appearance 
characteristic selector 14, an output media selector 16, an interpolation kernel 
lookup table memory 18, an interpolation processor 20, and a formatting 

10 processor 22. The system 10 further includes one or more output media provided 
by a cathode ray tube (CRT) display monitor 24 and a digital laser imager 26. The 
interpolation processor 20 receives image data, as indicated by block 28, from an 
input imaging device. In a medical imaging environment, the input imaging device 
may comprise, for example, a diagnostic device, such as a magnetic resonance 

1 5 (MR), computed tomography (CT), conventional radiography (X-ray), or 

ultrasound device. Thus, image data 28 may represent an image of a physiological 
object obtained for diagnostic purposes. 

The image data 28 received by interpolation processor 20 is representative 
of a plurality of pixels within the image. The number of pixels represented by 

20 image data 28 can be significantly less than the output resolution of output 
imaging media used for visible presentation of the image. The system 10 is 
configured to perform interpolation on image data 28 to produce interpolated 
image data representing the original pixels and interpolated pixels. The system 1 0 
selects one of a plurality of interpolation kernel lookup tables based on a 

25 combination of a selected format, a selected output media, and a selected 

appearance characteristic of the image. The interpolation kernel lookup tables 
contain interpolation coefficients that are applied to the original pixels in image 
data 28 to produce interpolated pixels, thereby transforming the spatial resolution 
of image data 28 to match the resolution of the output media. 

10 
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The format selector 12 selects one of a plurality of formats for 
presentation of an image on imaging media. The format selector 12 enables a 
system user to select a format representing a number of images to be formed on a 
angle page of the output imaging media. As an example, system 10 may provide 
5 formats for as many images as can be accommodated by a single page of media or, 
as a practical matter, as many images as can be viewed for useful diagnostic 
purposes. The term "page," as used herein, may refer to a sheet of photosensitive 
media such as film printed by laser imager 26, or a screen displayed by CRT 
display monitor 24. The format selector 12 can be realized, for example, by a user 
10 interface associated with a display panel on laser imager 26 or by a manual control 
panel on laser imager 26. The number of images specified by the selected format 
determines the size of each individual image to be formed on the imaging media. 
The size of each image determines the amount of expansion necessary to match 
the resolution of image data 28 to that of the output media. 

1 5 The appearance characteristic selector 14 selects one of a plurality of 

appearance characteristics of the presentation of the unage on the imaging media. 
The appearance characteristics may include a plurality of degrees of apparent 
sharpness of the image relative to a degree of apparent smoothness, when viewed 
by a human observer. The terms "sharpness" and "smoothness," as used herein, 

20 refer to degrees of apparent edge contrast between original pixels and interpolated 
pixels in the image, when viewed by a human observer. The appearance 
characteristic selector 14 enables a user to select a degree of apparent edge 
contrast in the visible representation of the image when formed on the imaging 
media to provide a more pleasing or diagnostically useful appearance. The format 

25 selector 12 can be realized, for example, by a user interface associated with a 

display panel on laser imager 26 or by a manual control panel on laser imager 26. 

The media selector 16 selects one of a plurality of imaging media for 
visible presentation of the image represented by image data 28. Specifically, 
media selector 16 enables the user to select other CRT monitor 24 or laser imager 

30 26, and also the type of imaging media printed by the laser imager. For example, 

11 
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laser imager 26 may be configured to print the image on either wet chemically 
processed silver halide film or thermographic dry silver film. The laser imager 26 
conceivably could be configured to print images on other types of thermographic 
film. The media selector 16 can be realized, for example, by a user interface 
S associated with a display panel on laser imager 26 or by a manual control panel on 
laser imager 26. Alternatively, media selector 16 may be configured to respond 
automatically to media selection information sent with image data 28 from an 
input imaging device. 

The interpolation kernel lookup table memory 18 stores a plurality of 

10 different interpolation kernel lookup tables. Each of the interpolation kernel 

lookup tables contains a plurality of interpolation coefficients calculated according 
to particular interpolation kernel functions. The interpolation kernel lookup tables 
may include, for example, at least one interpolation kernel lookup table containing 
a plurality of two-point interpolation coefficients and at least one interpolation 

15 kernel lookup table containing a plurality of four-point interpolation coefficients. 
In particular, memory 18 may store a variety of interpolation kernel lookup tables 
including an interpolation kernel lookup table containing a plurality of pixel 
replication interpolation coefficients, an interpolation kernel lookup table 
containing a plurality of bilinear interpolation coefficients, an interpolation kernel 

20 lookup table containing a plurality of third order polynomial interpolation 

coefficients, an interpolation kernel lookup table containing a plurality of fifth 
order polynomial interpolation coefficients, and one or more interpolation kernel 
lookup tables containing a plurality of cubic spline interpolation coefficients. 
The interpolation processor 20 selects one of the interpolation kernel 

25 lookup tables stored in memory 18 based on a combination of the format selected 
by format selector 12, the media selected by media selector 16, and the 
appearance characteristic selected by appearance characteristic selector 14. The 
particular interpolation kernel lookup table selected by interpolation processor 20 
produces interpolated image data having an optimum visible appearance, given the 

30 format, media, and appearance characteristic specified by the user. The 

12 
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interpolation processor 20 determines a sampling increment based on the amount 
of pixels required by the selected format and the amount of original pixels 
provided by image data 28. The interpolation processor 20 then applies 
interpolation coefficients contained in the selected interpolation kernel lookup 
5 table that correspond to the sampling increment to image data 28 to produce the 
interpolated image data. 

The format processor 22 receives the interpolated image data produced by 
interpolation processor 20, and prepares the interpolated image data for 
presentation on a single page of output imaging media. Specifically, format 

10 processor 22 assembles the interpolated image data produced for each of a 
number of images to be printed or displayed on a single page according to a 
format selected by the user. The format processor 22 assigns the interpolated 
image data for each image to a particular area of the page, and then rasterizes the 
assembled image data for application to drive either display monitor 24 or laser 

IS imager 26. 

Fig. 2 is a functional block diagram of the interpolation processor 20 
shown in Fig. 1, in accordance with the present invention. As shown in Fig. 2, 
processor 20 includes an input image data buffer 30, a row interpolator 32, a 
column interpolator 34, and an output image data buffer 36. The input image data 

20 buffer 30 receives image data 28 from an input imaging device. The row 

interpolator 32 processes original pixel values along each row of image data 28 
and interpolates them based on interpolation coefficients retrieved from the 
selected interpolation kernel table, as indicated by block 38. The row interpolator 
32 multiplies the original pixel values by corresponding interpolation coefficients 

25 38, selected according to the sampling rate, and sums the product to produce 
interpolated row pixel values. The resulting interpolated row pixel values each 
correspond to an interpolated pixel in the image, and represent an intensity of that 
interpolated pixel. The interpolated row pixel values create additional columns of 
interpolated pixels between columns of original pixels in image data 28. The 

30 intensity of each interpolated pixel is a function of the intensities of the original 

13 
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pixel values to which the interpolation coefficients were applied, and the values of 
the interpolation coefficients. 

Once a plurality of rows of original pixels have been interpolated by row 
interpolator 32, column interpolator 34 processes them to produce interpolated 
column pixel values. Specifically, column interpolator 34 interpolates both the 
original row pixels and the interpolated row pixels along a column based on 
interpolation coefficients 38 retrieved from the selected interpolation kernel table. 
The column interpolator 34 multiplies the original row pixels and the interpolated 
row pixels along a column by corresponding interpolation coefficients, and sums 
the product to produce interpolated column pixel values. The interpolated column 
pixel values create additional rows of interpolated pixels between rows of original 
pixels and interpolated pixels. Again, the intensity of each interpolated column 
pixel is a function of the intensities of the original column pixel values to which the 
interpolation coefficients were applied, and the values of the interpolation 
coefficients. After processing by row interpolator 32 and column interpolator 34, 
the interpolated image data is sent to output image data buffer 36. The image data 
is stored in output image data buffer 36 until retrieved by format processor 22. 

The row interpolator 32 and column interpolate 34 perform interpolation 

of image data 28 in a well known manner. The following FORTRAN code is 

provided, however, for purposes of illustration. The code illustrates the 

computation involved in a four-point interpolation process using a cubic spline 

interpolation kernel function, a two-point interpolation process using a fifth order 

polynomial function, and a two-point interpolation process using bilinear 

interpolation. 

c Four-point interpolation Cubic Spline 

c Parameter " ca " is set to " -1. " 

c ul, u2, u3, u4 are the calculated coefficients from 

c Cubic Spline kernel function 

c In this example, a 200% interpolation is shown. 

c 

ca = -l. 

do 100 H),4 

14 
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cl =float(iy4.0 

ul(i) - ca*(cl**3H2.*ca*(cl**2)>+ca*cl 
u2(i) = (ca+2)*(cl»*3) -(ca+3.)*(cl**2>* 1. 
u3(i) = -(ca+2.)*(cl**3)+<(2.*ci^3)*(cl**2))-ca*cl 
u4(i) = -ca*(cl**3>+«a*(cl **2) 
100 continue 



c 
c= 



c Two point interpolation using a 5th order polynomial 
10 c In this example, a 200% interpolation is shown 



c= 



do 200 i=0,4 
cl =float(i>/4.0 
ul(i) = 0. 

15 u2(i) =l-(10.*(cl**3))+(15.*(cl**4)H6.*(cl**5)) 

u3(i)-l-(10^(l-cl)**3))+(15.*((l-cl)»*4)H6.»((l-cl)**5)) 
u40) = 0. 

c 

c Two point interpolation using Bilinear 
20 c In this example a 200% interpolation is shown 



c= 



cl =float(i)/4.0 
ulCO^O. 
u2(i) = 1-cl 
25 u3C0 = cl 

u4(i) = 0. 
200 continue 



c= 



c In this example "mbb(i,k)" are the expanded image data 

30 c and "mnb(i,lc)" are the unexpanded origina] image data. 

c *************** interpolation p^w ************** 

c 

do 301 i = 1,4 
do 301 k = 0,3 

35 mbbftk+l)^I(k)*mnbai>+u2(k)*^^ 
1 +u4(k)*mnb0,4) 

if(mbb(i,k+l) .ge. 4095)mbb(i,k+l)=4095 
if(mbbO,k+l) le. 0)mbbO,k+l)=0 
301 continue 
40 c 

c***************Int«poIation Colurnn************* 
c 

do 302 i=l,4 
do302k=0,3 

45 mb(k+l,0^1(k)*mbb(l,i>+ii2(^ 
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1 4u4(k)*mbb(4,i) 

if(mb(k+U) .gc. 4095)mb(k+l,i)=4095 
if(mb(k*U) .le. 0)mb(k+l,i)=0 

302 continue 



c= 



An example of the interpolation kernel lookup tables stored in memory 1 8 
and the maimer in which interpolation processor 20 selects one of the interpolation 
kernel lookup tables will now be described. In this example, memory 18 stores 
interpolation lookup tables containing interpolation coefficients that represent the 

10 following interpolation kernel functions: 

A. Two-Point Interpolation Kernel Functions. 

1. Pixel Replication Kernel Function. This two-point interpolation 
kernel lookup table contains interpolation coefficients representative of a 
conventional pixel replication interpolation kernel function Ki(x), whereby 

IS interpolated pixels between two adjacent original pixels are simply assigned the 

intensity values of the closest original pixel. This interpolation kernel function can 
be expressed as: 



20 where x is a variable representing a location of the interpolated pixel relative to 
the original pixels. The pixel replication interpolation kernel produces a sharp 
apparent edge contrast between pixels in the interpolated image data. 

2. Bilinear Interpolation Kernel Function. This two-point 
interpolation kernel lookup table contains interpolation coefficients representative 

25 of a conventional bilinear interpolation kernel function K^x), whereby 

interpolated pixels between two adjacent original pixels are assigned linearly 
weighted intensity values based on distance from the closest adjacent original 
pixel. This interpolation kernel function K^x) can be expressed as: 



K,(x)=l 
Ki(x) = 0 



|x|<=0.5 
M>0.5 



30 



Kz(x)=l-|x| 
K2(x) = 0 



0 <= |x| <= 1 
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This bilinear interpolation kernel function produces a smoother apparent edge 
contrast between pixels in the interpolated image data, relative to pixel replication. 

3. Third Order Polynomial Kernel Function. This two-point 
interpolation kernel lookup table contains interpolation coefficients representative 
of a conventional third order polynomial kernel function Ka(x), whereby 
interpolated pixels between two adjacent original pixels are assigned intensity 
values based on the following expression: 

K 3 (x) = 2|xf - 3|xf +1 0 <= M <= 1 

K 3 (x) = 0, | X |>1 

This third-order polynomial kerne! function produces an apparent edge contrast 
between pixels in the interpolated image that is sharper than that produced by the 
bilinear interpolation kernel function K 2 (x), but smoother than that produced by 
the pixel replication interpolation kernel function Ki(x). 

4. Fifth Order Polynomial Kernel Function. This two-point 
interpolation kernel lookup table contains interpolation coefficients representative 
of a conventional fifth order polynomial kernel function K^x), whereby 
interpolation pixels between two adjacent original pixels are assigned intensity 
values based on the following expression: 

K4(x) = -6|x| 5 +15|x| 4 - 10|xp +1 0 < |x| < 1 

K4(x) = 0 | X |>1 
The fifth order polynomial kernel function produces an apparent edge contrast 
between pixels in the interpolated image that is sharper than that produced by both 
the third order polynomial interpolation kernel function K 3 (x) and the bilinear 
replication interpolation kernel function K 2 (x), but smoother than that produced 
by the pixel replication interpolation kernel function K t (x). 
B. Four-Point Interpolation Kernel Functions. 

5. Cubic Spline Function (a = 0). This four-point interpolation kernel 
lookup table contains interpolation coefficients representative of a conventional 
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cubic spline kernel function Kj(x), whereby interpolation pixels between two 
adjacent original pixels are assigned intensity values based on the following 

expression: 

K 5 (x) = (a + 2)|x| 3 - (a + 3)|xf +1 0 <= |x| <- 1 
5 K 5 (x)=a|x| 3 -5a|x| 2 + 8a|x|-4a |x|>l 

K 5 (x) = 0 |x|>2 
With the coeffident a = 0, the above interpolation kernel function Ks(x) provides a 
relatively smooth apparent edge contrast between pixels in the interpolated image. 
6. Cubic Spline Function (a = -.2). This four-point interpolation 
10 kernel lookup table corresponds to the above expression Ks(x) with the coefficient 
a - -0.2 to realize an interpolation kernel function K*(x). With a = -0.2, the 
interpolation kernel function K*(x) produces an apparent edge contrast between 
pixels in the interpolated image that is sharper than that produced by the 
interpolation kernel function K 5 (x). 
IS 7. Cubic Spline Function (a = -.4). This four-point interpolation 

kernel lookup table corresponds to the above expression K 5 (x) with the coeffident 
a = -0.4 to realize an interpolation kernel function Kt(x). With a = -0.4, the 
interpolation kernel function K^x) produces an apparent edge contrast between 
pixels in the interpolated image that is sharper than that produced by the 
20 interpolation kernel function K«(x). 

8. Cubic Spline Function (a = -.6). This four-point interpolation 
kernel lookup table corresponds to the above expression K 5 (x) with the coeffident 
a - -0.6 to realize an interpolation kernel function K*(x). With a = -0.6, the 
interpolation kernel function Ks(x) produces an apparent edge contrast between 

25 pixels in the interpolated image that is sharper than that produced by the 
interpolation kernel function Kt(x). 

9. Cubic Spline Function (a = -.S). This four-point interpolation 
kernel lookup table corresponds to the above expression K 5 (x) with the coeffident 
a - -0.8 to realize an interpolation kernel function K^x). With a - -0.8, the 

30 interpolation kernel function K*(x) produces an apparent edge contrast between 
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pixels in the interpolated image that is sharper than that produced by the 
interpolation kernel function Kt(x). 

10. Cubic Spline Function (a = -1 .0). This four-point interpolation 
kernel lookup table corresponds to the above expression K s (x) with the coefficient 

5 a = -1 .0 to realize an interpolation kernel function Kio(x). With a = -1.0, the 
interpolation kernel function Kto(x) produces an apparent edge contrast between 
pixels in the interpolated image that is sharper than that produced by the 
interpolation kernel function Kg(x). 

1 1 . Cubic Spline Function (a = -1.2). This four-point interpolation 

10 kernel lookup table corresponds to the above expression K 5 (x) with the coefficient 
a = -1.2 to realize an interpolation kernel function K n (x). With a = -1.2, the 
interpolation kernel function Kn(x) produces an apparent edge contrast between 
pixels in the interpolated image that is sharper than that produced by the 
interpolation kernel function Ki 0 (x). 

15 12. Cubic Spline Function (a = -1.4). This four-point interpolation 

kernel lookup table corresponds to the above expression K 5 (x) with the coefficient 
a = -1 .4 to realize an interpolation kernel function K 12 (x). With a = -1 .4, the 
interpolation kernel function K t2 (x) produces an apparent edge contrast between 
pixels in the interpolated image that is sharper than that produced by the 

20 interpolation kernel function K u (x). 

13. Cubic Spline Function (a = -1 .6). This four-point interpolation 
kernel lookup table corresponds to the above expression K 5 (x) with the coefficient 
a = -1.6 to realize an interpolation kernel function K| 3 (x). With a = -1 .6, the 
interpolation kernel function K»(x) produces an apparent edge contrast between 

25 pixels in the interpolated image that is sharper than that produced by the 
interpolation kernel function K i2 (x). 

14. Cubic Spline Function (a = -i .8). This four-point interpolation 
kernel lookup table corresponds the above expression K 3 (x) with the coefficient a 
= -1.8 to realize an interpolation kernel function K u (x). With a = -l .8, the 

30 interpolation kernel function K I4 (x) produces an apparent edge contrast between 
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pixels in the interpolated image that is sharper than that produced by the 
interpolation kernel function K»(x). 

15. Cubic Spline Function (a = -2.0). This four-point interpolation 
kernel lookup table corresponds to the above expression K 5 (x) with the coefficient 

5 a = -2.0 to realize an interpolation kernel function K !5 (x). With a = -2.0, the 
interpolation kernel function Kis(x) produces an apparent edge contrast between 
pixels in the interpolated image that is sharper than that produced by the 
interpolation kernel function Km(x). 

1 6. Cubic Spline Function (a = -2.2). This four-point interpolation 

10 kernel lookup table corresponds to the above expression K 3 (x) with the coefficient 
a = -2.2 to realize an interpolation kernel function Kj 6 (x). With a = -2.2, the 
interpolation kernel function K ]6 (x) produces an apparent edge contrast between 
pixels in the interpolated image that is sharper than that produced by the 
interpolation kernel function K»(x). 

1 5 The interpolation processor 20 selects one of the above interpolation 

kernel lookup tables Ki-Km to improve the appearance of a visible representation 
of the interpolated image data based on the selected format, the selected media, 
and the selected appearance characteristic. The selection is made by interpolation 
processor 20 in recognition that the results of an interpolation kernel function can 

20 vary depending on the degree of expansion or reduction required by the format 

selected by format selector 12, and the response characteristic of the output media 
selected by output media selector 14. The processor 20 also bases the selection 
on an appearance characteristic, e.g., sharpness or smoothness, selected by 
appearance characteristic selector 14. In this manner, interpolation processor 20 

25 selects from memory 1 8 the best available interpolation kernel function Ki-K I6 for 
the present conditions. 

In this example, interpolation processor 20 first narrows the field of 
interpolation kernel lookup tables stored in memory 18 according to the format 
selected by format selector 12. For selected formats requiring smaller expansion 

30 rates, interpolation processor 20 selects one of the two-point interpolation kernel 
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lookup tables K1-K4, in view of the lower incidence of blurring with small-scale 
expansion. Because the smaller expansion rate produces less blurring, a two-point 
interpolation kernel lookup table can provide acceptable visible results. In this 
manner, the selection of a two-point, rather than four-point, interpolation kernel 
5 lookup table by interpolation processor 20 reduces the amount of computational 
effort necessary to produce the interpolated image data. The computation time 
required for two-point interpolation will be less than that required for four-point 
interpolation by at least a factor of four. For formats requiring higher expansion 
rates, interpolation processor 20 selects one of the four-point interpolation kernel 

10 tables K4-K16. The interpolation processor 20 thereby produces better spatial 
frequency when needed to avoid the higher incidence of blurring that can occur 
with larger expansion. The four-point interpolation kernel requires added 
computation, but is necessary to ensure a pleasing or diagnostically useful visible 
presentation of the interpolated image data on the output media. 

IS As an example of the criteria applied by interpolation processor 20 to 

select a type of interpolation kernel lookup table, satisfactory visible results have 
been observed for the use of two-point interpolation kernel tables with selected 
formats requiring expansion rates that do not exceed a factor of three in both 
horizontal and vertical directions. For selected formats that require expansion 

20 rates of greater than a factor of three in both horizontal and vertical directions, 
interpolation processor 20 should be configured to select one of the four-point 
interpolation kernel tables. At expansion rates of less than a factor of three, the 
degree of blurring that may occur with a two-point interpolation kernel has been 
found to be acceptable, particularly in view of the computational savings gained. 

25 At expansion rates of greater than a factor of three, the incidence ofblurring 
generally dictates the use of an interpolation kernel producing better spatial 
frequency. Thus, the improved visible results achieved by use of a four-point 
interpolation kernel are considered to be worth the added computational effort. 
Of course, the question of whether the viable presentation of the interpolated 

30 image data is acceptable will ultimately be a subjective matter for the human 
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observer. Therefore, other criteria for selecting between a two-point or four-point 
interpolation kernel are conceivable. 

After interpolation processor 20 has narrowed the selection process to 
either the set of two-point interpolation kernel lookup tables or the set of four- 
5 point interpolation kernel lookup tables, the process turns to selection of a subset 
of lookup tables appropriate for the type of output imaging media selected by 
output media selector 16. Both the set of two-point interpolation lookup tables 
K1-K4 and the set of four-point interpolation kernel lookup tables K5-K16 include 
interpolation coefficients representing interpolation kernel functions designed to 

1 0 produce different appearance characteristics, i.e., different degrees of sharpness or 
smoothness. For example, the set of two-point interpolation kernel tables 
includes, in descending order of apparent sharpness, the pixel replication 
interpolation kernel table K u the fifth order polynomial interpolation kernel table 
the third order polynomial interpolation kernel table K 3 , and the bilinear 

1 5 interpolation kernel table K 2 . Similarly, the set of four-point interpolation kernel 
tables includes, in descending order of apparent sharpness, cubic spline 
interpolation kernel tables K 5 -Ki 6 modified by the value of coefficient a to adjust 
sharpness. 

The interpolation processor 20 selects a subset of appropriate interpolation 
20 kernel lookup tables by considering the response characteristics of the output 
imaging media, which can greatly affect the visible results of the interpolation. 
Thermographic dry silver film has a development speed that is significantly slower 
than that of conventional silver halide film, in both toe and shoulder speeds. The 
scanning laser in laser imager 26 may have an insufficient dynamic range to 
25 compensate for the slower development speed. Consequently, blurring can result 
in high contrast areas of the viable presentation of the interpolated image data, 
such as those areas involving a transition from image data representing 
alphanumeric text at a minimum density to image data representing pictorial 
information at much higher densities approaching maximum density. Similarly, a 
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CRT display monitor suffers from severe dynamic range limitations, and therefore 
can also exhibit blurring in high contrast areas. 

When the media selected by output media selector 16 is conventional silver 
halide film, interpolation processor 20 selects one of the interpolation kernel 
5 lookup tables K r K 16 producing a relatively smooth apparent edge contrast. The 
interpolation processor 20 selects an interpolation kernel lookup table producing a 
relatively smooth apparent edge contrast in recognition of the fester development 
speed of the silver halide film, and hence the smaller dynamic range required to 
reproduce the desired optical densities. As an example, if a selected format 

10 requires an expansion rate of less than three in both horizontal and vertical 

directions, and the selected media is silver halide film, interpolation processor 20 
selects one of two-point interpolation kernel tables K1-K4 providing a relatively 
smooth apparent edge contrast, such as the bilinear interpolation kernel table K 2 . 
In contrast, if a selected format requires an expansion rate of greater than three in 

1 5 both horizontal and vertical directions, and the selected media is silver halide film, 
interpolation processor 20 selects one of the four-point interpolation kernel tables 
K3-K16 providing a relatively smooth apparent edge contrast. For example, the 
sixth four-point interpolation kernel table K*, with coefficient a set to -0.2, 
provides a particularly smooth apparent edge contrast when applied to the 

20 interpolated image data. 

When the media selected by output media selector 16 is a thermographic 
dry silver film, interpolation processor 20 selects one of the interpolation kernel 
lookup tables Ki-Km producing a sharper apparent edge contrast than that for 
silver halide film. The interpolation processor 20 selects an interpolation kernel 

25 lookup table producing a relatively sharp appearance in recognition of the slower 
development speed of the dry silver film, and hence the larger dynamic range 
required to reproduce the desired optical densities. The sharper interpolation 
kernel lookup table provides a mechanism for improving the visible appearance of 
the presentation of the interpolated image data on the dry silver film. Specifically, 

30 the sharpness produced by the selected interpolation kernel lookup table produces 
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an "overshoot/undershoot" effect that compensates for the slower development 
speed of the dry silver film. However, the overshoot and undershoot are clipped 
in high contrast areas by the minimum density Dmin and the maximum density 
Dmax. 

The sharpness of the selected interpolating kernel lookup table increases 
the apparent spatial frequency of the interpolated image data by increasing the 
apparent slope between adjacent interpolated pixels. In particular, the selected 
interpolating kernel lookup table increases the apparent slope between minimum 
density Dmin and maximum density Dmax in high contrast areas, such as those 
involving a transition from image data representing alphanumeric text at a 
minimum density to image data representing pictorial information at much higher 
densities approaching maximum density. Of course, the other areas of the visible 
presentation having varying densities also will be affected by the sharper 
interpolation kernel lookup table. The effect is not readily apparent, however, for 
interpolation with formats requiring smaller expansion factors. Further, the four- 
point interpolation kernel lookup tables selected for formats requiring larger 
expansion fectors conceal the effect by avoiding the introduction of a large 
amount of high spatial frequency artifact and aliasing. 

As an example, if a selected format requires an expansion rate of less than 
three in both horizontal and vertical directions, and the selected media is dry silver 
film, interpolation processor 20 selects one of the two-point interpolation kernel 
tables K1-K4 providing a relatively sharp apparent edge contrast, such as the pixel 
replication interpolation kernel table Ki, the third order polynomial kernel table 
K 3 , or the fifth order polynomial kernel table K4. In contrast, if a selected format 
requires an expansion rate of greater than three in both horizontal and vertical 
directions, and the selected media is dry silver film, interpolation processor 20 
selects one of four-point interpolation kernel tables K 5 -Ki 6 providing relatively 
sharper apparent edge contrasts. For example, the interpolation kernel table Ku, 
with coefficient a set to -1 .8, provides a very high degree of overshoot and 
undershoot, producing a very sharp apparent edge contrast when applied to the 
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interpolated image data. In areas in which pixels transition from Dmin to Dmax, 
however, the overshoot and undershoot are clipped. Thus, the sharper 
interpolation kernel lookup table provides the benefits of higher apparent slope 
between Dmin and Dmax, compensating for the slow development speed of the 
5 dry silver film, but does not suffer significantly from the overshoot and undershoot 
effect, which would ordinarily be undesirable. 

When the media selected by output media selector 16 is CRT display 
monitor 24, interpolation processor 20 selects one of the interpolation kernel 
lookup tables K r Ki 6 producing an even sharper apparent edge contrast than that 

10 required by the dry silver film. The interpolation processor 20 selects an 

interpolation kernel lookup table producing a relatively sharp appearance in 
recognition of the severely limited dynamic range of the CRT monitor. As with 
the dry silver film, the sharper interpolation kernel lookup table provides a 
mechanism for improving the visible appearance of the presentation of the 

1 5 interpolated image data on the CRT monitor. Specifically, the sharpness produced 
by the selected interpolation kernel lookup table produces the overshoot and 
undershoot effect to provide a higher apparent slope between maximum intensity 
and minimum intensity on the monitor. The overshoot/undershoot effect 
compensates for the limited dynamic range of the CRT monitor, but is clipped in 

20 high contrast areas by the minimum and maximum intensities. The sharpness of 
the selected interpolating kernel lookup table increases the apparent spatial 
frequency of the interpolated image data by increasing the apparent slope between 
adjacent interpolated pixels in high contrast areas, such as those involving a 
transition from image data representing alphanumeric text at a minimum intensity 

25 to image data representing pictorial information at much higher intensities 
approaching maximum intensity. 

As an example, if a selected format requires an expansion rate of less than 
three in both horizontal and vertical directions, and the selected media is CRT 
monitor 24, interpolation processor 20 selects one of the two-point interpolation 

30 kernel tables K1-K4 providing a sharp apparent edge contrast. The pixel 
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replication interpolation kernel table Ki and the fifth order polynomial kernel table 
K4, in particular, produce interpolated image data having sharper apparent edge 
contrasts. If a selected format requires an expansion rate of greater than three in 
both horizontal and vertical directions, however, and the selected media is CRT 
5 monitor 24, interpolation processor 20 selects one of the four-point interpolation 
kernel tables K5-K16 providing relatively sharper apparent edge contrasts. For 
example, the interpolation kernel table K t6 , with coefficient a set to -2.2, provides 
a maximum degree of overshoot and undershoot, producing a very sharp apparent 
edge contrast what applied to the interpolated image data. The sharper 
10 interpolation kernel lookup table provides a higher apparent slope to compensate 
for the limited dynamic range of CRT monitor 24, thereby improving the visible 
presentation of the interpolated image data to a human viewer. 

After interpolation processor 20 has narrowed the selection process to 
either the set of two-point interpolation kernel lookup tables K r Ki or the set of 
1 5 four-point interpolation kernel lookup tables K 5 -K 16 based on the selected format, 
and to a particular subset of the interpolation kernel lookup tables Ki-Ki 6 based 
on the selected media, the process turns to selection of a particular interpolation 
kernel lookup table based on the appearance characteristic selected by appearance 
characteristic selector 14. Thus, even though the selection by interpolation 
20 processor 20 is largely driven by format and media considerations, it still may be 
desirable to allow a system user some latitude for sharpness and smoothness. As a 
first example, if a selected format requires an expansion rate of greater than three 
in both horizontal and vertical directions, and the selected media is dry silver film, 
interpolation processor 20 selects one of the four-point interpolation kernel tables 
25 K3-K16 providing a sharp apparent edge contrast. The apparent sharpness 

produced by each of interpolation kernel tables K 5 -Ki 6 increases in response to an 
increase in the coefficient a. Thus, if h is determined that interpolation kernel 
table K< produces a sufficient overshoot/undershoot effect to improve the visible 
presentation of the interpolated image data, but the system user still desires a 
30 sharper appearance, interpolation processor 20 selects a higher-numbered 
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interpolation kernel lookup table having a coefficient a that corresponds to the 
desired sharpness. 

Fig. 3 is a flow diagram illustrating a method for adaptive interpolation of 
image data, in accordance with the present invention. It will be recognized that 
the method of the present invention can be readily implemented by the system 10 
described above. As shown in the flow diagram of Fig. 3, the method includes a 
step of selecting one of a plurality of formats for presentation of an image on an 
imaging media, as indicated by block 40. A particular imaging media is selected 
for presentation of an image, as indicated by block 42. Further, one of a plurality 
of appearance characteristics, representing sharpness or smoothness, are selected 
for the presentation of the image on the imaging media, as indicated by block 44. 
As indicated by block 46, one of a plurality of different interpolation kernel lookup 
tables then is selected based on a combination of the selected format, the selected 
media, and the selected appearance characteristic. Each of interpolation kernel 
lookup tables contains a plurality of interpolation coefficients representing the 
interpolation carried by a different interpolation kernel functions, as described 
above with respect to system 10. Finally, as indicated by block 48, the 
interpolation coefficients contained in the selected interpolation kernel lookup 
table are applied to the image data to produce the interpolated image data. 

EXAMPLE 

The following example is provided to illustrate a system and method for 
adaptive interpolation of image data, in accordance with the present invention. In 
particular, the following example illustrates one possible implementation of the 
selection process carried out by the system and method of the present invention. 

In this example, one of interpolation kernel lookup tables Kj-K* is selected 
if the selected format requires an expansion factor of less than or equal to three in 
horizontal and vertical directions, and one of interpolation kemei lookup tables 
Ks-Kn is selected if the selected format requires an expansion factor of greater 
than three in horizontal and vertical direction. 
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If the selected media is silver halide and the selected fbnnat requires an 
expansion factor of less than or equal to three, one of interpolation kernel lookup 
tables K1-K4 is selected according to the selected appearance characteristic, with 
Ki representing the sharpest appearance and K 2 the smoothest. If the selected 
5 media is silver halide and the selected format requires an expansion factor of 

greater than three, one of interpolation kernel lookup tables in the subset K 5 -K a is 
selected according to the selected appearance characteristic, with K 5 representing 
the smoothest appearance in the subset and K X2 the sharpest. 

If the selected media is dry silver and the selected format requires an 

10 expansion factor of less than or equal to three, one of interpolation kernel lookup 
tables K1-K4 is selected according to the selected appearance characteristic, with 
Ki representing the sharpest appearance and K 2 the smoothest. If the selected 
media is dry silver and the selected format requires an expansion factor of greater 
than three, one of interpolation kernel lookup tables in the subset K*-K u is 

15 selected according to the selected appearance characteristic, with K* representing 
the least sharp appearance in the subset and K u the sharpest. 

If the selected media is CRT monitor 24 and the selected format requires 
an expansion factor of less than or equal to three, one of interpolation kernel 
lookup tables K1-K1 is selected according to the selected appearance 

20 characteristic, with Ki representing the sharpest appearance and K 2 the smoothest. 
If the selected media is CRT monitor 24 and the selected format requires an 
expansion factor of greater than three, one of interpolation kernel lookup tables in 
the subset Ki 0 -Ki 6 is selected according to the selected appearance characteristic, 
with Kio representing the least sharp appearance in the subset and K 16 the 

25 sharpest. 

Having described the exemplary embodiments of the invention, additional 
advantages and modifications will readily occur to those skilled in the art from 
consideration of the specification and practice of the invention disclosed herein. 
Therefore, the specification and examples should be considered exemplary only, 
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with the true scope and spirit of the invention being indicated by the following 
claims. 
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What is claimed is: 

1. A method for performing interpolation on image data 
representative of a plurality of pixels within an image to produce interpolated 
image data, said method comprising the steps of: 

5 selecting one of a plurality of formats for presentation of said image on 

imaging media; 

selecting one of a plurality of different interpolation kernel lookup tables 
based on the selected format, each of said interpolation kernel lookup tables 
containing a plurality of interpolation coefficients; and 
10 interpolation means for applying the selected interpolation kernel lookup 

table to said image data to produce said interpolated image data. 

2. The method of claim 1 , wherein said plurality of interpolation 
kernel lookup tables includes at least one interpolation kernel lookup table 
containing a plurality of two-point interpolation coefficients and at least one 

1 5 interpolation kernel lookup table containing a plurality of four-point interpolation 
coefficients. 

3. The method of claim 1, wherein said plurality of interpolation 
kernel lookup tables includes an interpolation kernel lookup table containing a 
plurality of pixel replication interpolation coefficients, an interpolation kernel 

20 lookup table containing a plurality of bilinear interpolation coefficients, an 

interpolation kernel lookup table containing a plurality of third order polynomial 
interpolation coefficients, an interpolation kernel lookup table containing a 
plurality of fifth order polynomial interpolation coefficients, and an interpolation 
kernel lookup table containing a plurality of cubic spline interpolation coefficients. 

25 4. The method of claim 1, further comprising the step selecting one of 

a plurality of appearance characteristics of said presentation of said image on said 
imaging media, said step of selecting one of said interpolation kernel lookup tables 
further including a step of selecting one of said interpolation kernel lookup tables 
based on both the selected format and the selected appearance characteristic. 
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5. The method of claim 4, wherein said plurality of appearance 
characteristics includes a plurality of degrees of apparent sharpness of said image 
relative to a degree of apparent smoothness of said image. 

6. The method of claim 1, further comprising the step of selecting one 
5 of a plurality of imaging media for presentation of said image, said step selecting 

one of said interpolation kernel lookup tables including selecting one of said 
interpolation kernel lookup tables based on both the selected and the selected 
media. 

7. The method of claim 6, wherein said plurality of media includes a 
10 silver halide film, a thermographic film, and a cathode ray tube monitor. 

8. The method of claim 1, wherein said image is one of a plurality of 
images, each of said plurality of formats specifying a number of said images on 
said imaging media, said number of said images determining a size of each of said 
images. 

IS 9. A method for performing interpolation on image data 

representative of a plurality of pixels within an image to produce interpolated 
image data, said method comprising the steps of: 

selecting one of a plurality of imaging media for presentation of said 

image; 

20 selecting one of a plurality of different interpolation kernel lookup tables 

based on the selected media, each of said interpolation kernel lookup tables 
containing a plurality of interpolation coefficients; and 

interpolation means for applying the selected interpolation kernel lookup 
table to said image data to produce said interpolated image data. 

25 1 0. The method of claim 9, wherein said plurality of interpolation 

kernel lookup tables includes at least one interpolation kernel lookup table 
containing a plurality of two-point interpolation coefficients and at least one 
interpolation kernel lookup table containing a plurality of four-point interpolation 
coefficients. 
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11. The method of claim 9, wherein said plurality of interpolation 
kernel lookup tables includes an interpolation kernel lookup table containing a 
plurality of pixel replication interpolation coefficients, an interpolation kernel 
lookup table containing a plurality ofbilinear interpolation coefficients, an 

5 interpolation kernel lookup table containing a plurality of third order polynomial 
interpolation coefficients, an interpolation kernel lookup table containing a 
plurality of fifth order polynomial interpolation coefficients, and an interpolation 
kernel lookup table containing a plurality of cubic spline interpolation coefficients. 

12. The method of daim 9, further comprising the step selecting one of 
10 a plurality of appearance characteristics of said presentation of said image on said 

imaging media, said step of selecting one of said interpolation kernel lookup tables 
including selecting one of said interpolation kernel lookup tables based on both the 
selected media and the selected appearance characteristic. 

13. The method of claim 12, wherein said plurality of appearance 

1 5 characteristics includes a plurality of degrees of apparent sharpness of said image 
relative to a degree of apparent smoothness of said image. 

14. The method of claim 9, wherein said plurality of media includes a 
silver halide film, a thermographic film, and a cathode ray tube monitor. 
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